GLOBAL SMALL SOLUTIONS TO THE CRITICAL 
RADIAL DIRAC EQUATION WITH POTENTIAL 



FEDERICO CACCIAFESTA 



Abstract. We solve globally a radial cubic Dirac equation perturbed with a 
small potential, with data of small critical norm H^. The main tool are new 
endpoint estimates of the perturbed Dirac flow for a class of radial-type initial 
data. 

1. Introduction 
Consider a zero mass nonlinear Dirac equation 

iut-Vu = F{u), u{0,x)^f{x) (1.1) 
for the spinor field m : Mj x Ki^ ^ C*, where V is the operator defined by 

3 
k=l 

while the 4x4 Dirac matrices are defined as 

a.= (^0 f^-W (1.2) 

in terms of the Pauli matrices 

'^i = ( J J ) , = ( ° 7 ) , as = J \ y (1.3) 

We shall assume that the nonlinear term F{u) is cubic of a very specific form, 
namely 

either F(u) — {/3u,u)u or F(u) — {u,u)u, (1-4) 
where (•, •) denotes the standard hermitian product in and f3 is the 4x4 Dirac 
matrix 

■ I2 



-I. 

Notice that nonlinearities of the form (1.4) are the most interesting from a physical 
point of view (see e.g. [15]). 

The typical approach to the Dirac equation is based on the identity 

{idt - v){idt + v) = i-dt + A)i4 (1.5) 

which follows from the anticommuting relations 

atak + akcti = 2Sikl4, i,fc = l,2,3, (1.6) 
a,/3 + /3ai = 0, i = 1,2,3, ^^=14. (1.7) 
In many cases the identity allows to reduce problems for the massless Dirac equation 
to analogous ones for the wave equation, for which many effective tools are available. 



Date: January 13, 2013. 

2000 Mathematics Subject Classification. 35J10, 35Qxx, 42B20, 42B35 . 

Key words and phrases, singular integrals, weighted spaces, Schrodinger operator, Schrodinger 
equation, Strichartz estimates, smoothing estimates. 

1 



CRITICAL DIRAC EQUATION 



2 



The classical strategy to solve nonlinear dispersive equations is by the use of a 
fixed point argument in a suitable space, via the appropriate space-time Strichartz 
estimates. A huge literature is available on these estimates for several dispersive 
operators (wave, Schroedinger, Klein-Gordon, Dirac and others). For homogeneous 
nonlinear terms one typically finds a threshold regularity Sc in the scale of Sobolev 
spaces ^ such that for subcritical data with s > Sc solvability holds, while for 
supercritical data with s < Sc one has various degrees of ill-posedness. Data of crit- 
ical regularity give rise to difficult questions which depend on the precise structure 
of the equation. 

If we denote with e''^^ the propagator for the Dirac operator V, which can be 
defined via Fourier transform as the operator of symbol e**^ ", Strichartz estimates 
for the linear 3D Dirac equation can be written as follows (see [9], [10], and, for the 
Dirac equation, [2], [.3]): 

\\\D\^-^^-h^'''fh.L.<\\fh2 (1.8) 

where the exponents (p, q) are admissible, that is to say 
2 2 

- + - = 1, 2<p<oo, 2<q<oo. 
P q 

Here and in the following we shall use the notation \D\ = (— A)^/^, the mixed space- 
time Strichartz space L| is LP{M.t; £''(]R")), and is the homogeneous Sobolev 
space with norm = || |I?|*/||l2. As it is well known, estimate (1.8) is false 

for the endpoint couple {p,q) = (2,oo), which would correspond to the estimate 

l!e"^/llL?L,=o< 11/11^,. (1.9) 

The failure of (1.9) for generic data was first noticed for 3D the wave equation in 
[11], and the corresponding statement for the Dirac equation follows easily from 
the relation connecting the Dirac and wave flows (see (2.2) below). 

A clever use of the estimates (1.8) is sufficient for the study of (1.1) in the 
subcritical case. Local existence, and global existence for small data, was proved 
in [7] for nonlinearities of the form (/3w, with p > 3 and data in with 

s > 3/2 — 1/ip— 1). The case of cubic nonlinearities and small data in , s > 1, 
was solved in [13], where the result was generalized to all space dimensions. We 
also mention that systems involving the Dirac equation, like Dirac-Klein-Gordon 
and Maxwell-Dirac, have beed the object of intense attention and the subcritical 
theory was completed only recently (see [5], [4], [(>]). 

The failure of the endpoint estimate (1.8) means that above methods break down 
in the critical case of a cubic nonlinearity with data, and indeed the problem of 
global existence is still open in this case. In an attempt to overcome this limitation, 
in [12] the following precised estimate was proved: 

V''^f\\LlL^^^Ll<\\f\\m Vp<(X). (1.10) 

The norm at the left hand side distinguishes between the integrability in the radial 
and tangential directions. Using estimate (1.10), Machihara et al. were able to 
prove global well posedness for problem (1.1) for small iJ^-norm data with slight 
additional angular regularity, and in particular for all radial data. This is 
especially interesting since, as we shall see, radial data do not correspond to radial 
solution for the Dirac equation. 

The main goal of the present paper is to study a special class of solutions to 
equation (1.1) and more generally to a potential perturbation of the form 



iut-Vu + V{x)u^ F{u), u{0) = f{x), 



(1.11) 



CRITICAL DIRAC EQUATION 



3 



which can be considered as a suitable generahzation of the radial solutions to the 
wave equation in the context of the Dirac equation. We recall that, for small 
potentials V with suitable decay at infinity, the full range of Strichartz estimates 
(1.8) holds also for the perturbed fiow e'*^'^''"^^ as proved in [-3]. Hence subcritical 
problems can be treated exactly as in the unperturbed case and one can extend the 
results of [7], [13] to (1.11) in a straightforward way. Here we focus on the more 
difficult case of critical data with an additional symmetry assumption like 

/ - /i + 2?/2 with hEH\ h(.H\ /i, /2 radial. 

In addition, in order to preserve the symmetry of solutions, we need to assume that 
the potential V{x) is spherically symmetric in the sense of [IG]. 
Our first result is an endpoint estimate for the linear flow: 

Theorem 1.1. Let V{x) he a Ax A matrix of the form 

V{x) ^Vi{\x\)li+il3{a-x)V2{\x\), Vi,V2 ^ M., x G (1.12) 

(where x — x/\x\). Assume that for some a > \ and some sufficiently small 5 > Q 

Then the following endpoint Strichartz estimate 

l|e''^''+'''/llL?Ls= < WfWm (1-14) 
holds for all initial data f E where 

V} = {/i +P./2, /i e ij^(M3), /2 e ij2(M3), /i,/2 radial}. (1.15) 

We recall that condition (1.13) is sufficient to ensure that the perturbed Dirac 
operator V + V is self-adjoint on the domain _ff^(R3)^. For this and many other 
properties of the Dirac equation, a comprehensive reference is [IG]. 

The natural application of estimate (1.14) is to prove global well posedness for 
the critical equation (1.11). However the nonlinear term F(u) does not operate on 
the space if} and additional restrictions on the algebraic structure of the data are 
necessary. More precisely, it is possible to decompose the space L'^{E?Y a direct 
sum 

oo j 

L2(r3)4^ L2(o^+oo;dr)®iJ„^,fc^.. 

2 2 ±(3 + 1/2) 

where the spaces Hmj,kj are two dimensional and are generated by spherical har- 
monics on the sphere (this is called a dcomposition in partial wave subspaces). 
When j — 1/2, we have four spaces 

L^{0, +oo; dr) (S) H^^^^^k,,2 

corresponding to the four possible choices of indices 

(TOi/2,fci/2)- (-1/2,-1), (-1/2,1), (1/2,-1), (1/2,1). (1.16) 

Then we notice the mildly surprising fact that each of these four spaces is invariant 
not only for the Dirac operator but also for the action of cubic nonlinearities of the 
forms (1.4). In a sense, these spaces can be considered as a suitable generalization of 
radial functions adapted to the structure of the nonlinear problem (1.11). A detailed 
analysis of the partial wave decomposition is given in Section 5, with explicit forms 
for the functions in these spaces (see Lemma 5.5 and in particular (5.22)-(5.25)). 
Thanks to this invariance, we can prove the following global existence result: 
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Theorem 1.2. Consider the equation on M x M'^ 

iut-Vu + V{x)u^ F{u), u{0,x)^f{x) (1.17) 
where the potential has the form 

V = V^i(|x|)l4 + ^/3(a • x)V2{\x\), x = x/\x\ 

and satisfies assumption (1.13), while 

F{u) ~ {l3u,u)u or F{u) = {u,u)u. 

Assume the initial data f belong to a space H^{{0,oo),dr) ®^mi/2.fei/2 ""^^ '^f 
the choices (1.16). Then if the norm of the data is sufficiently small, problem 
(1.17) has a unique global solution in the class Ct{M.,H^) f] Lf {R,L°°). 

The paper is organized as foUows. Sections 2 contains an extension of the end- 
point estimate for the free Dirac operator which is then adapted in Section 3 to 
a mixed endpoint-smoothing estimate with weights for the nonhomogeneous hnear 
Dirac equation. Section 4 is devoted to the proof of Theorem (1.1). In Section 5 
we recall the structure of the Dirac operator, the partial wave decomposition, and 
we investigate the interaction of the algebraic structure with the nonlinear term. 
Section 6 contains the proof of Theorem (1.2). 

The author wishes to thank prof. Piero D'Ancona for his help and the referee 
for some useful suggestions. 



2. The homogeneous endpoint estimate 

Consider the 3-dimensional, massless, free Dirac equation (1.1). As already 
observed, the Dirac flow does not preserve radiality so that we cannot hope to 
repeat here the simple argument used in [11] for the wave equation. However we 
can prove an endpoint estimate for suitable classes of function; to this end we need a 
deeper insight in the structure of the Dirac operator expressed in radial coordinates. 

Let u = e**'^/ and notice that, thanks to identity (1.5), u is (formally) a solution 

of 

'□m = 

uiO) = fix) (2.1) 
[ut{0) = iVf. 



This gives the formula 



:^*^/ = cos(t|C|)/ + .^^^^P/ (2.2) 



and we easily see that this representation is valid for generic distribution data. 
We start by proving the following result: 

Proposition 2.1. Let f belong to the space T-L^ defined in (1.15). Then the follow- 
ing endpoint Strichartz estimate holds: 

l|e'*^/llL?L,== < II/IIh- (2-3) 

By formula (2.2), we see that the proof is an immediate consequence of the following 
Lemma. Notice that the proof of the first estimate (2.4) is inspired by an argument 
of [8] which holds for all n > 3 without modification, while we fix n = 3 in the 
second estimate. 
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Lemma 2.2. Let f be a radial function and let e'*'^' be the linear propagator 
associated to the wave operator. Then the following estimates hold: 

l|e'*l^l/liL?Ls= < 11/11^^ (">3) (2.4) 

t\D\ 



D ^ 



<\\f\\m (« = 3). (2.5) 



Proof. Using Fourier transform in spherical coordinates and the radiahty of / and 
setting p = |a;|, = A, a; • ^ = pAcos6', we have 



OO /'TV 



^^t\D\J^l^{x■^+m)j^^^^^^^ / / e*^(*+''™"''V'(A)A"-i(sin6i)"-2rf6'dA. 
J Jo Jo 

(2.6) 

With the change of variable y = cos 6* in (2.6) we obtain 

dAe"^.g(A) J e'^^^ (1 - y^)"^ dy 

with g{X) = /(A)A"~^i?(A) {H represents the classical Heaviside function). Now 
changing the order of the integrals we obtain 

/ + 00 />1 
e^Pit+^y)g(p)dp= / dy{l-y')'^g{t + Xy) 
-OO J -I 

Since for n > 3 one has (1 — y^)~^ < 1, the change of variable y ^ yjr yields 
again 

<- r g{t + y)dy^M{g)it) (2.7) 

^ J —r 

where M denotes the standard maxhnal operator. Then we have for all t 

l|e^*l''l/llL^<M(g)(t), (2.8) 

and thus by the L^-boundedness of maximal operator and Plancherel's theorem 

1 1 

l|e^*l''l/llL?L,=o <||5Lj= ^ (X-'lflfdXj =^ A^l/I A-^dAj = 

(2.9) 

= ||A"^/|U. = 11/11^:^ 

which gives (2.4). 

We now turn to estimate (2.5), for which the calculations are similar. Indeed we 
can write (here we are fixing n — 3) 

^Vf = J e'(-€+*l€l)(a • 0/(0^^ = (2.10) 

where 

X sr^ («fc ■ Cfc) 

St 

Using spherical coordinates as before we have 

dX d<l) dee'^^'+P''°'''^'>A{e,(l))f{X)X^ sine (2.11) 
/o Jo Jo 

with the operator A{6, (f)) = ai cos + a2 sin 9 cos 4> + sin 9 sin 0. Observing that 
a2 sin 9 cos (f) d(f) — / as sin sin c?(/) = 0, 
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we see that (2.11) is equal to 



dX / rf6le'^(*+''™"'')aiCos6l/(A)A2sin6i. 
Jo 

Setting as before g{X) — /(A)A^iJ(A), changing variable cos 9 y and then y y/r 
yield 

dy {a, -y) dp e^^(*+''^)g(p) = - / (a^■^ g{t + y) dy = cM{g){t) 

■1 J-oo T J~r ^ ' 

(2.12) 



since the term ^ai • — ^ is bounded, and so we have the bound 



ait\D\ 



-Vf 



< 



M{g){t)- 



(2.13) 



The LP-boundcdness of maximal operator and Placherel Theorem yield as above 

it\D\ 

< WfWm (2.14) 



-Vf 



D 



which gives the desired estimate (2.5). 



□ 



Combining estimates (2.4) and (2.5) and using representation (2.2) for the solu- 
tion of the free Dirac system we obtain estimate (2.3). 

3. The mixed endpoint-smoothing estimate 
We consider now the non homogeneous equation 

iut-Vu^ F{t,x), u{Q,x)^0. (3.1) 
By Duhamel's formula and the representation (2.2) we can write the solution u as 



u{t,x) = / e*(*-'')^F(s,x)ds 



(3.2) 





* ^cos((i-s)|D|)F(s,2;) +i^^^ifc-f^^2?F(s,x) ) ds. 



\D\ 

Thus in order to estimate the solution u to (3.1) we can deal separately with the 
two integrals 

e'(*-^)'^'f^(s,x)ds and / ^— 



-VF{s,x)ds. 



We prove the following: 

Proposition 3.1. Let n — i and assume F{t,x) has the structure 

F{t,x) ^ Fi{\x\)l4 + ip{a ■ x)F2{\x\). (3.3) 
Then the following estimate holds 



f e<t--)T^F{s) ds 
Jo 



< \\(x)^+\D\F\ 



LiLl 



(3.4) 



The key step in the proof of (3.4) is the following non homogeneous estimate for 
the wave propagator with a radial term. 

Lemma 3.2. Let n > 3, F{t,-) be a radial function. Then the following estimate 
holds 



< 



{xr^ + \D\^FU2tLl. 



(3.5) 
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Proof. We start with (3.5). Expanding u as in the homogeneous case (see formulas 
(2.6 and (2.7)), from the radiality of F we can estimate the L°° norm of the solution 
at fixed t as (here G{s, A) — A"~^F(s, X)H{X) and H is the Heaviside function) 



^ sup 



Gis,y + t-s) 



ds\ dy = 



supiy (^J G{s,y + t~s){y + t-s)^+{y + t-s)^'^'dsjdy 



< 



sup ■ 



< 



dy 



Gi{s,y + 1- s) 



2 \ 2 

ds 



{y + t~ sy^^ds 



where in the last inequality we have used Cauchy-Schwarz inequality and Gi is the 
function defined by 

Gi{s,y) = G{s,y){y)^+. 



ds we have 



Setting now h{z) :— [ j Gi{s,z — s) 
1 



sup ■ 



dy 



Giis,y + t-s) ^ ds]' = M{h){t) 



The boundedness of the maximal operator yields 



Gi{s,t-s) 



dsdt 



\Gi\ 



The last quantity is precisely 

\\{y)i+Tx^y[r-'F{s,X)H{X)) 

and to conclude the proof we need to estimate it by 

<\\{x)i + \D\^FU.^.. 

Since T~^{\D\^~ f) = \^\^~ f we see that it is enough to prove the general in- 
equality (we can neglect the dependence on time) 

(A^/(A)i/(A)) 11^. < ||(x)'^7IIl^ (3.6) 

for k = 1/2+ . 

We prove (3.6) by interpolation. The case A: = is trivial, since from Placherel's 
Theorem we obviously have 

fA^/(A)i7(A)) - ||A^/(A)|U2 = ll/IU.. 



(3.7) 

The case fc = 1 is just a little more complicated. Since of course (p) < 1 + \p\, we 
need only prove that 



(3.8) 



\\pT,^, [X^f{X)HiX)) |1l^<||(x)/|U. 

or equivalently 

\\d^[x^f{X)H{X)) <||(x)/|U.. 
We write 

ii^A (a^7(a)xe+(a)) iu^ < wx^dxfimLi + u'^fWhi = h+h 

with the shorthand notation L\ = L^{Q, oo). For Ii we trivially have 

Il = \\X^(^)\\L2 9i\\\x\f\\L2. 
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Let's now turn to I2. We split the norm 

l|A^/(A)L^ = \\\^f{\)\\m{x>i}) + ||A^/(A)||i.({A<i}). 
Plancherel's Theorem yields again for the first term 

I|A'^/(A)||l2({a>i}) < I|A'^/||l2(r) 11/11^2. 
To handle the second term we use Hardy's inequality: 

l|A'^/(A)|U.(^,<i}) = |||erV1lLMiei<i}) ^ ^ IIN/IU^ 

and this concludes the case k = 1. By interpolation with (3.7) we obtain the desired 
estimate (3.5). □ 

Lemma 3.3. Let n — i and F(t,-) be of the form (3.3). Then the following 
estimate holds 

ft ^i(t-s)\D\ 



\D\ 



-VFis) ds 



< 



{x)^+\D\Fh2tLl- 



(3.9) 



for k = 0,1. 



Proof. Since the operator i'Dj3{a ■ x)(j) applied to a radial function (j) produces the 
radial function iPcf)" , Lemma (3.2) holds, and we only need to control terms of the 
form (where F^ad denotes a radial function) 



t ^i(t-s)\D\ 



AjFrad{s) ds 



with Aj = i[3^{a ■ x), j = 0, 1. Recalling (2.10)-(2.12), since the quantities Aj are 
obviously bounded, we can estimate in both cases with 



t g»(t-s)|D| 



-AjFradis)ds 



< 



* 1 



\G{s,y + t-s)\dsdy 



where as before G{s, A) = X^Frad{s, X)H{X). Proceeding exactly as in the proof of 



Lemma (3.2) we obtain estimate (3.9). 



□ 



Estimate (3.4) is an immediate consequence of (3.5), (3.9) and representation 
(3.2) 



4. Proof of Theorem (1.1) 

We now turn to the proof of Theorem (1.1). This is based on a simple application 
of a smoothing estimate for a Dirac equation with potential proved in [3] (see also 
[2], [1] for related results): 

Theorem 4.1 ([3]). Let V{x) — V{x)* be a 4x4 complex valued matrix and assume 
that for some cr > 1 and some sufficiently small 6 > one has 

S 



\Vix)\ < 



w„{x) 

Then the following smoothing estimate holds 



where w„{x) = \x\{l + | log |a;| 



iw^^'-e-- ' ■ VllLji^ < ||/||l2. 



(4.1) 



(4.2) 



It is not difficult to deduce the endpoint estimate (1.14) for the perturbed flow 
from the previous result and our mixed endpoint-smoothing estimate (3.4). First 
of all, the solution of the equation 

iut ~ Vu + Vu, w(0, x) = f 
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can be written, regarding Vu as a right-hand member of the equation 

Jo 

Then we can write 

\\\D\-^u\\^.r^ = |l|i?|-ie^*(^+^)/Lj,. < (4.3) 



^-^(T/(.)e-(^+^)/)d. 



The first term can be estimated by (2.3) (notice that \D\ commutes with V and 
hence with the flow). In order to apply estimate (3.4) to the second term we need 
the following 

Lemma 4.2. If f E , then e**'^/ G , and if V is of the form (1.12), then 
Ve'^'^f is of the form (3.3). 

Proof. We write f = fi +2^/2 with /i, /2 radial functions. Then we have, from 
(2.2), 

e-/ = (cos(.li^l) + ^^I^) if^+Vf,) = 
= cos{t\D\)f, + sm{t\D\)\D\f2+v(cos{t\D\)f2 + ""fflf^'V i 



\D\ 

where /i and /2 are radial functions with the appropirate regularity, and this con- 
cludes the proof of the first statement. The proof of the second statement is triv- 
ial. □ 

We thus can estimate (4.3) with (2.3) and (3.4) obtaining 

<II/IIl^ + II(^)^+^«IIl?l^ 

Now multiplying and dividing by w^ixY^^ in the second norm on the right hand 
side yields 

< ii/iu. + \\{xy/^+wi/'vh^ ■ \\w-'/\\mLi- 

Notice that the weighted norm of V at the right hand side is bounded as it fol- 
lows from assumption (1.13). Moreover (1.13) implies also that the assumption of 
Theorem 4.1 is satisfied. Then using (4.2) we conclude 

\\\D\-\h.^^ < (l + ||(a;)V2+^yV|U.) Wfh. (4.4) 

that gives, under hypothesis (1.13) on the potential, the desired Strichartz endpoint 
estimate. 



5. Partial wave subspaces and radial Dirac operator 

The purpose of this section is to construct, following [l(i], invariant subspaces 
for the Dirac operator with a potential having a special symmetry. To this end 
we use the classical decomposition of the space L^(M'^)'' in the direct sum of 2- 
dimensional Hilbert spaces, the partial wave subspaces, which are invariant for the 
Dirac operator. We shall aslo check that the lowest order partial wave subspaces 
are invariant even for the cubic nonlinearities that we consider here. 

We begin by recalling the basic facts on the decomposition, referring to [IG] for 
more details. We shall use the standard notation for polar coordinates in R"^ 

X ^ r sin 6 cos 0, y = r sin 6 sin cf), z — r cos 9 
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with the unit vectors in the directions of the polar coordinate hnes given by 

X 

Cr = (sin 9 cos 4>, sin 6 sin 0, cos 9) = — : = x 
eg — (cos 9 cos (j>, COS 9 sin <j), — sin 9) = — ^ 

1 T 

60 = (-sin(/i,cos0,O) = 

Then we write for a function V' G L'^{^^) 

^{r, 9, 0) = ri^{x{r, 9, c^),v{r, 9, </)), z(r, 9, 0)). (5.1) 

Since the function ^p{r,-,-) of the angular variables is square integrable on the unit 
sphere L^(S'^), the mapping ip ^ ip defines a unitary isomorphism 

L2(M3) ^ L^{{0, oo), dr; L''{S^)) = L\{0, oo), dr) (g, L''{S^). 

Applying the transformation (5.1) on each component of the (vector valued) wave- 
function, we obtain the analogous decomposition 

i2(M3)4 ^ L^(^(^o^oo),dr) (g, L^iS^. 

The decomposition of the Hilbert space into a " radial" and an " angular" part is 
very useful since the angular momentum operators 

Jj = X A (— iV) orbital angular momentum 
J — L + S total angular momentum 
act only on the angular part L?{S'^)'^ in a nontrivial way; here 

S = -l/4(aAa) 

is the spin angular momentum operator. Recalling the expression of V in polar 
coordinates 

d I ( d I d\ 

V = e,— + - 667^ + 60 — — 5.2 
or r \ o9 sm 9 o9 I 



we obtain that, since x = 



r ■ Cr 



L = zee^— — - 160— 5.3 
sm 9 a<p o9 

where the differentiation applies to each component of the wavefunction. 

The Dirac operator can be written in polar coordinates as follows. Combining 
(5.2) and (5.3) yields 

. d I, 

—iV = —I e^-r (6r A L) 

or r 

and thus 

d 1 

— ia • V = —iia ■ Cr)-:: a ■ icr A L). (5.4) 

Or r 

By the basic property of the Dirac matrices: 

[a ■ A){a ■ B) = A- B + 2iS- (AAB); 

which holds for any matrix-valued vector fields A — {A^^A^^A^), B — {B^,B'^, B^) 
with F\G' e 7W4x4(C), and 

75Q; = 2S, 

where 75 — ^ ]' Q ) ' obtain 

(a • A){2S ■ B) = ^75^1 ■ B-ia-{AhB), 
thus equation (5.4) is equal to 

d i 

= —i{a ■ Cr)— + -(a • 6r)(2S • L). 
or r 
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Finally, introducing the spin orbit operator 

if = /3(2S-L + 1) = ^(J2-L2 + i/4) (5,5) 

where we used the identity — {L + S)"^ — + 2S ■ L + 3/4, we arrive at the 
following representation: 

Proposition 5.1. The 3- dimensional Dirac operator can be written as 

V = -i{a-x)(-^ + - --fSK] (5.6) 
\or r r J 

where K is the spin orbit operator defined in (5.5). 

The key step to construct the invariant spaces is the following: 

Proposition 5.2. For each choice {j,mj,kj) with j — nij = + 

kj — — (j + 1/2), +(j + 1/2), there exist precisely two eigenfunctions 
k- ^ C°°{S'^)^ satisfying the following relations: 

K^mj^kj = -kj^mj,kj- 

The family <i>^ j, forms an orthonomral basis of L^[Ei^)^. 

The functions ^mj.kj can be written explicitly using spherical harmonics. We 
first recall the following representation of 3-diinensional spherical harmonics 



'(^''^^ = \/^^^|tW^""'^"('^°^^) V - / < m < L (5.7) 
where P™ are the Legendre polynomials 

As it is well known, the spherical harmonics form a complete orthonormal set in 
L^(5^), i.e. every function / G L^(5^) can be written as 

00 / 

for some constants /™; moreover, they are eigenfunctions of both the operators 
and L3, i.e. 

i2y-™ ^ 1(^1 + i^Yi"' (5.9) 
LsYi"" = mYi"'. (5.10) 
We now define for j — i, |, |, .., mj = — j, — j + 1, + j the functions ^'^'^1/2 ^ 
L2(52)2: 

*-V2-^|^^/J^^™.|l/2 j (5-11) 



"^^■+^2 - V2JT2 I, -V^- + i + „,,y™^jV2 j ■ (5.12) 

These functions are, as it is easily seen, eigenfunctions of both the operators 
and J2 = L'^ + a ■ L + 3/4 with eigenvalues l{l + l) and j(j + 1) respectively. So we 
conclude, in view of (5.5), that the functions in Proposition (5.2) are given by 

/ i*"' \ / \ 
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Thus the Hilbert space L^(S'^)'* is the orthogonal direct sum of 2-dimensional 
Hilbert spaces 'Hmj,kj, which are spanned by simultaneous eigenfunctions ^. 
ofJ^andX: 

oo j 

L'iSr= (5-14) 

j=^,l,...m,=-jk,=±U+^} 

Easy calculations show that the functions 5'J2i/2 satisfy 

(a ■ x)^""' = 



and hence 



*(«-*X..,. =T$^w,fc,- (5.15) 



This proves the following: 



Lemma 5.3. The subspaces 'Hmj,kj o-i"e left invariant by the operators P and 

a ■ X. With respect to the basis kj^mj kj} defined above, these operators are 

represented by the 2x2 matrices 

The decomposition just shown obviously implies a similar one of L^(R'^)'*, in 
which each partial wave subspace L^{{0, oo), dr)®'H„ij,kj is isomorphic to i^((0, oo), drY 
if we choose the basis {^^. /.^ , j,^ }• There is in fact a unitary isomorphism be- 
tween the Hilbert spaces: 

2,2(^3)4 ^ ^2((o^ ^ ^^^^^^ 

This decomposition and (5.6) allow us to easily calculate the action of the Dirac 
operator (at least on differentiable states) even in the presence of a suitable poten- 
tial. 

Proposition 5.4. The Dirac operator (5.6) with the potential 

V{x) = Vi{\x\)li^^ip{a-x)V2{\x\) (5.18) 

leaves the partial wave subspaces C^(0, oo) ^Hmj.kj invariant. With respect to the 
basis <i>^. fc. the Dirac operator on each subspace can be represented by the 

operator 

d r^(l^l) -^r + '-F + VM) \ .5 19) 

"-'^"U + v + ^^dxD v,{\x\) ) ^'-''^ 

which is well defined over C|5"(0,oo)2 c L^((0,oo),dr)2. Moreover, the Dirac oper- 
ator D on (M.'^)^is unitary equivalent to the direct sum of the partial wave Dirac 
operators djnj,kj, 

oo j 

rf™.,-^, (5-20) 
i=i,|,...™j=-jfc,=±o+i) 

Remark 5.1. Proposition 5.4 holds for slightly more general potentials (see [16]), 
but we shall not need this fact here. 

Remark 5.2. The operator in (5.19) is also known as radial Dirac operator. It can 
be proved that dmj,kj is essentially self-adjoint (for every j) on C^(0, oo) if and 
only if P + is essentially self-adjoint on C(f'(R3\{0}). 
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Thus using spherical coordinates it is possible to construct invariant spaces for 
the perturbed Dirac operator. What may come as a surprise is that for j = 1/2 the 
partial wave subspaces are also invariant for the cubic nonlinearity, and this fact is 
obviously crucial for the nonlinear application we shall prove in the next section. 

Lemma 5.5. Let j = 1/2 and let (^1/2,^1/2) one of the couples (-1/2,-1), 
(-1/2,1), (1/2,-1), (1/2,1). Then the partial wave subspaces {{0, 00), dr) (g) 
^mi/2,fci/2 '^''^ invariant for the cubic nonlinearities Fi{u) — {u,u)u and F2{u) = 
u)u, i.e. 

C(f ((0,oo),dr)®H™,^,,fc,,, ^Fi{u), F2{u) G C^{{Q,^),dr) ®n^^^^^k,;,- 

(5.21) 

Proof. We explicitly write down the functions $+, $^ in the four cases: a straight- 
forward calculation using formulas (5.7), (5.8), (5.11), (5.12) and (5.13) yields 



-1/2,-1 



2^ 


V ) 
( 



-1/2,-1 



20F 



e'-^sinf 



cos B 



1/2,1 



20F 



2V^ 





1/2,1 



V 20F / 



(5.22) 



(5.23) 



1/2,-1 



1/2,1 



/ 



20F 



V / 



1/2,-1 



1 



2^ 

V ^e*'^sin6' 
\ 2^ 



( 



: COS 



2V^ 

* e*-* sin 61 



20F 



V 



1/2,1 



* 

1 

V / 



(5.24) 



(5.25) 



We prove Lemma (5.5) for the couple (1/2, 1), i.e. for functions of the form (5.25), 
being the proof for the other cases completely analogous. 

The generic function u E i^((0, 00), dr) (E) ^1/2,1 can be written as 

uir, 9, 0) = u+ ir)^+/2,ii^, cp) + (r)*-/^^, [9, 0) 

for some radial functions , . So / takes the vectorial form 



( u'^{r) „ ^ ^ cos 9 \ 



2^ 



[r) 



2^ 
u (r) 

2^ 




e"^ sin( 



(5.26) 
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Thus the Hcrmitian product (u, u) yields 

{u,u) = — — cos^ 6 u+{rf — sin^ 9 u+{r f + —u^{rf = 

47r 47r An 

that has no angular eomponcnts. This proves that if u G C^((0, oo), dr)<S)'Hmi/2,ki/-y 
then g Co°°((0,oo),dr) 

Minor modifications yield the same result also for the nonlinearity F2 (u) . In fact 
we know from Lemma (5.21) that the operator /3 acts on the partial wave subspaces 
in a very simple way with respect to the basis 'I'^}: if in fact we associate to 
the function u its coordinates (u^(r), M~(r)) with respect to such a basis we have 
/3u = (u"'"(r), — u~(r)), so that 

{I3u,u) = -^cos^ eu+{rf - ^sin^ eu+{rf - -^^"(r)^ = 
47r 4tt 4tt 

that again has no angular components, and this shows that if u e C^{{0, 00), dr) (E) 
^mi/2,fei/2 thenF2(u) € C^iiO,<x),dr) (g,Hrn,^,.M/2- ^ 

6. Global existence for the nonlinear equation 

As an application of the results we have presented in the previous section we can 
now prove global existence for problem (1.17) with small initial data in one of the 
four partial wave subspaces H^{{0, 00), dr) 'Hmi/2,fci/2- Ou'" go^-l to prove 

Theorem 6.1. Consider the Cauchy problem for the 3-dimensional nonlinear Dirac 
equation 

iut-Vu + V{x)u^ F{u), u{0,x) = f{x) (6.1) 

where the potential V is of the form 

V = Vi{\x\)l4 + i/3{a ■ x)V2{\x\) 

and satisfies assumption (1.13 ), while the nonlinear term F(u) is either of the form 
{l3u, u)u or {u, u)u. 

Then for every initial data f G H^{{0, 00), dr) ® 'Hmi/^M/^' with sufficiently 
small norm, there exists a unique global solution u(t,x) to problem (1-17) in 
the class Ct{R,H^)nL^t{R,L°°). 

Proof. The proof is identical for both choices of the form of the nonlinear term. We 
rewrite (1-17) in integral form 

u = e**^/ + i f e*(*-")^ {V{s)u{s) + F{u{s))) ds = (6.2) 
Jo 

= e"^/ + « f e''^'-'^^ iV{s)u{s))ds + i f e'<^'^'^^ {F{u{s))) ds ^ 
Jo Jo 

= e**(^+^)/ + i f e'(*-^)^ {F{u{s))) ds = $(u) = h + h 
Jo 

we denote by $(u) the RHS of (6.2) and we check that the map <i> is a contraction 
on the function space 

X^LlL^f^LrHl. 
In order to estimate the first term Ii we use our endpoint Strichartz estimate 
(1.14), observing that if / G i?^((0, 00), dr) ®T~Lmi/2M/2 then in particular / is of 
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the form f = fi + 'T)f2 with /i, /2 radial functions, so estimate (1.14) holds and 
gives 

l|/i||x< 11/11^. (6.3) 
Now we need to handle the nonlinear term 12- By Minkowski inequality 



e*(*-'')^F(u(s))ds 



X 



By energy conservation we have 

On the other hand, in view of Lemma (5.5), we can use estimate (2.3) and we have 

||e^*^e-^F(u(.))L.^. = ||^^(^(.))||^. < \\F{u{s))\\^^^,. 

Thus 

\\h\\x<\\F{u{s))\\^^^, 
Then by Holder inequality in t, x we obtain 



u\ 



which implies 

ll<i'(")IU<ll/lli^i + lkllx- 

An analogous computation gives 

mu)-^v)\\x<{\\u\\l + \\v\\l)\\u-v\\x. (6.4) 

Therefore if the data belong to a sufficiently small ball in H^, $ is a contraction 
on that ball, and its unique fixed point is the unique global solutions to problem 
(1.17) in the space X. □ 
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